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We report results of our calculation on the i:Na term and quark spin content of the nucleon on the quenched 
16^ X 24 lattice at /3 = 6.0. The disconnected insertions which involve contributions from the sea quarks are 
calculated with the stochastic Z'z noise algorithm. As a physical test of the algorithm, we show that the forward 
matrix elements of the vector and pseudoscalar currents for the disconnected insertions are indeed consistent with 
the known results of zero. Wc tried the Wuppertal smeared source and found it to be more noisy than the point 
source. With unrcnormalizcd m,j = 4.42(17)MeV, we find the nNa term to be 39.2 ± 5.2 MeV. The strange 
quark condensate in the nucleon is large, i.e. (7V|ss|7V) = 1.16 ± 0.54. For the quark spin content, we find A-u = 
0.78 ± 0.07, Ad = -0.42 ± 0.07, and As = -0.13 ± 0.06. The fiavor-singlet axial charge = AE = 0.22 ± 0.09. 



1. Substantial sea effects 




The Recent SMC experiment measured gi 
structure function in polarized deep inelastic 
muon-proton scattoring which confirmed a large 
negative strange quark polarization in its contri- 
bution to the proton's spin [1] 
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-0.12 ± 0.04 ± 0.04 
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An analysis of tt N scattering also suggests a 
large strange quark contribution [2] 

{N I ss I N) 



, >~0-2 (2) 
{N \uu-\-dd + ss\ N) ^ ' 

Both of these results suggest that bilinear 
strange quark operator may have large matrix 
elements between proton states. In view of the 
fact that the strange is not a valence quark of the 
proton, both these results indicate that the 'sea' 
quark effects are substantial. 

The study of -kNo term also involves the sea 
quark effect and has been controversial over the 
years. Recent chiral perturbation calculation 
which takes into account the form factor effect in 
the extrapolation from the Cheng-Dashen point 
to where the -kNo term is defined gives a smaller 
value [3] 



G = mq{N I MM + I iV) ~ 45MeV 



(3) 



In principle, lattice QCD should be a direct 
way of accessing these sea quark effects through a 
numerical simulation. That is what we have done 
in the last several years [4,5] and we will report 



Figure 1. Sea quark loop in nucleon 



some of our results. We should point out that we 
missed a factor of 4 in our previous calculation of 
sea quark matrix element of the scalar currents 
[5] With this factor of 4, our preliminary results 
[5] are consistent with the final numbers to be 
reported here. 

2. Sea quark matrix elements on the lattice 

The major computational difficulty is calculat- 
ing the summation of quark loops over all space 
points at each time slice in order to obtain the 
zero momentum transfer (forward) matrix ele- 
ments. To this end, the stochastic Z2 noise al- 
gorithm is employed [6]. 

To extract the matrix element, we calculate the 
ratio of three- to two-point functions 
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Hence, the matrix element can be obtained as 
the slope from the ratio. 

3. Physical test of the algorithm 

Since this is the first application of the noise al- 
gorithm to the three point functions for the study 
of tho soa quark effects, it is essential to be able 
to show that the algorithm passes some tests and 
that the algorithm gives the correct physical sig- 
nals. To this end, we apply the Zi noise algorithm 
to calculate the disconnected quark loops of the 
vector and pseudoscalar form factors. It turns 
out that these form factors are proportional to 
the three momentum q. 

(A^ I r„ • VY"" I N) CX e^kradk 

(iV I r„ ■ P'°°^ \N)^qm 

Thus their forward matrix elements are zero. 

We carry out calculations for quenched Wilson 
fermion for 6 currents (local and point-split vector 
and axial, scalar, and pseudoscalar) on a 16^ x 24 
lattice with (3 = 6.0. Four quark masses with hop- 
ping parameters k = 0.154, 0.152, 0.148, 0.140, 
are studied. 300 stochastic steps are used for 
3 Hghter quarks and 200 steps are used for the 
heavier quarks. The calculations are done for 
24 gauge configurations, each separated by 1000 
Monte Carlo sweeps. 

As an example, the results for k = 0.148 are 
shown in Fig. 2. By examining the slopes of these 
curves, it is clear that the forward matrix ele- 
ments for the vector and pseudoscalar currents 
arc indeed consistent with zero. Given the suc- 
cessful test of the algorithm, we can proceed to 
estimate the wNa term and quark spin content of 
nucleon with more confidence. 

4. Wuppertal smearing technique 

We also tried the Wuppertal smearing tech- 
nique to improve our 2-point and 3-point function 
calculations. The Jacobi method is used to solve 
the equations and we take [8] Kg = 0.181, A'' = 
100 which corresponds to r ~ 4.1 for the 24 con- 
figurations. In our calculation with the Wup- 
pertal source, the signals appear earlier than the 
point source as expected. However, the error bars 
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Figure 3. The scalar results with Wuppertal 
source k = 0.148 



become larger. This can be seen in Fig. 3 as com- 
pared to the same case with the point source as 
shown in the upper- right corner of Fig. 2. It has 
been pointed out by the UKQCD group [8] that 
the Wuppertal smearing technique increases the 
fluctuation between configurations. We simply 
demonstrated this in the case of the disconnected 
insertion. 

5. Conclusion and discussion 

The minimum method is employed to fit the 
slopes of the ratios which are the matrix elements. 
Keeping the sea quark mass in the strange quark 
range (i.e. at k = .154) and extrapolating the va- 
lence quark mass to the chiral Hmit, we thus ob- 
tain the strange quark sea matrix elements. We 
used the fixed boundary condition in time, hence 
the t we sum ranges from the nucleon source at 
t = and sink at i = 16 to avoid the bound- 
ary contamination. Since the lowest state, i.e. 
the nucleon, appears at f > 8, the linear curves 
are fitted with t > 8. The mean-field improved 
one-loop renormalization is employed to renor- 
malize the matrix elements [9] . In the case of the 
the strange sea quark matrix elements, we need 
to consider the finite ma correction in addition. 
This has been reported by Lagae in this confer- 
ence [10]. We use the correction factors found in 
his study. The renormahzed results are reported 
as follows: 

<T = 39.2(5.2) MeV(unrenormalized mq = 4.42(17)MeV), 
{N\ss\N) = 1.16 ±0.54, 
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^oint iiplit vec:or currjent- 
M=-0.016(29) 
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Figure 2. Sea quark loop matrix elements within nucleon with k = 0.148 
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(iV I ss I N) 



{N \ UU + dd + ss \ N) 

2{N \ss\N) 



0.14 ±0.05, 



, , =0.33 ±0.09, 

{N \uu + dd\N) 

Audis = Ad,i,s = -0.13 ± 0.07, 

Alt = 0.78(7), Ad = -0.42(7), As = -0.13(6), 

pi = AE = 0.22(9), pi = 0.61(10), pi = 1.20(10). 

The TrNa term is consistent with the TriV scat- 
tering calculation [3] as well as the recent calcu- 
lation [7] at p = 5.7. The negative As is quite in 
agreement with the latest SMC experiment and, 
together with sea contributions from the u and d 
quarks, are large enough to cancel out the con- 
nected insertion part of the flavor-singlet axial 
charge to give a small g\ to be in agreement with 
the experimental finding. We believe this is a ma- 
jor step in resolving the "proton spin crisis" . We 
should emphasis that our calculation is gauge in- 
variant, hence the axial current does not fix with 
the gauge non-invariant dimension 3 gluon oper- 
ator. 
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